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ABSTRACT 

We consider the evolution of the probability distribution P(x, x, t), associated with 
an inhomogeneous light scalar field x i n the Robertson- Walker Universe, where the 
inhomogeneities are produced by quantum fluctuations during an earlier inflationary 
epoch. For a specific choice of scalar potential which occurs in models of so called 
late-time phase transitions in which domain walls are produced, P is shown to evolve 
from a Gaussian to a non-Gaussian distribution. The structure of the latter justifies 
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the recent use of 3-dimensional percolation theory to describe the initial distribution 
of domain walls in these models. 
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l. Introduction 



The idea that large scale structure in the Universe might appear as the result of a 
'late-time' phase transition associated with an extremely light scalar field, has been the 
subject of a number of recent papers [1,2]. Large scale structure in these models could 
be generated either by density perturbations due to fluctuations in this scalar field, 
or because of the gravitational attraction of topological defects such as domain walls 
that might be produced. If the walls form a network, then in order to avoid generating 
too large a microwave background anisotropy the phase transition should occur after 
photon decoupling. (Alternatively, the phase transition could occur prior to this time 
and still be consistent with observations of the microwave background if, for example, 
the domain walls that form are isolated 'bags' that subtend small angular distances 
on the sky.) For the phase transition to occur after photon decoupling requires the 
mass of the scalar to be extremely small, of order 1CT 29 eV, and for it to have very 
weak interactions [1]. Although masses of this size might at first sight appear to be 
unnatural, they have been shown to occur in several particle physics scenarios [2,3]. 

It should be understood that by 'phase transition' we do not necessarily mean 
one which occurs at finite temperature, which would normally require the scalar field 
responsible for wall formation to be in thermal equilibrium. Since we are dealing with 
a weakly interacting field, this latter requirement might be problematic (but see ref. [1] 
for a further discussion of this point.) In fact, the fluctuations in the scalar field that 
eventually produce the walls need not be thermal in origin, but might be the result of 
quantum fluctuations during an earlier inflationary epoch [4] . This is the viewpoint we 
shall adopt in this paper. 

It has been known for some time [5,6], that significant clustering of galaxies occurs 
at length scales of order IMpc. Moreover, detailed redshift surveys [7], have revealed 
the presence of walls, voids and other gravitational attractors on scales up to 50Mpcs, 
and at redshifts z < 5 which corresponds to relatively recent times in the Universe's 
evolution. It is possible that such clustering might be explained by primordial fluc- 
tuations alone, such as those produced during an earlier inflationary epoch. (In fact, 
such primordial perturbations are at least consistent with the recent discovery of a mi- 
crowave background anisotropy [8], <5T/T ps 6 x 1CT 6 ). However, it is difficult for such 
fluctuations to simultaneously account for structure on both large and small scales. A 
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better approach might be a kind of hybrid model containing primordial density per- 
turbations plus those due to the formation of late-time domain walls, where the latter 
would describe the kind of structures discovered by the redshift survey [7]. 

In refs.[3,9], a detailed analysis of the distribution of domain walls in such a model 
was presented, where the walls are formed by quantum fluctuations in the scalar field 
during an earlier inflationary epoch. Even though these fluctuations arise from in- 
flation, they do not generate walls until much later on, when inflation ends and the 
Robertson- Walker (R-W) phase begins. (The idea that domain walls can be produced 
in this way was first investigated, in a different context, in ref.[4]). An important as- 
sumption of the analysis in refs.[3,9], was that 3-dimensional percolation theory [10], 
provides a good description of the domain wall network immediately after their forma- 
tion. The criterion for percolation theory to be strictly applicable, is that the scalar 
field x describing the walls must be in either of two degenerate minima in all causally 
connected regions of space, with probability p and (1-p) respectively. By contrast, 
the probability distribution P(x,x) °f X immediately after inflation ends is a Gaus- 
sian centred about the value x- Since P(x,X) is the distribution of possible x values 
throughout space, it is clear that we have to do more to meet the above criterion. 

The purpose of this letter is to justify the use of percolation theory by showing that 
the distribution P(x, X) evolves during the R-W phase to a non-Gaussian distribution 
whose maxima approach two degenerate minima of the scalar potential. In this way 
the new distribution for x comes much closer to satisfying the percolation criterion 
discussed above. 

The paper is organized as follows : in section 2 we briefly review some concepts 
concerning scalar field quantization in de Sitter space and the appearance of the dis- 
tribution P(x, x). In section 3, we show how P(x,x) evolves into a non-Gaussian 
distribution and illustrate this with various numerical plots. 

2. Scalar field fluctuations in de Sitter space 

In this section we will briefly review some aspects of scalar field quantization in de 
Sitter space, and explain how this leads to the production of an inhomogeneous quasi- 
classical field. The reader is referred to refs.fll] for details. First consider quantization 
of a scalar field in the de Sitter geometry ds 2 = dt 2 — a 2 (t)dx 2 , where x and t are 
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( x > = / 7r^72 { °p e ip ' x P (t) + h.c } (2) 



comoving coordinates. The scale factor a(t) = e Hit , and H; is Hubble's constant during 
inflation. Q The equation of motion for is 

4> - a- 2 (t) V V + 3fL> + ^ = (1) 

d0 

Note that we choose the time variable t = t — t; where t is actual cosmological time 
lapsed, and t; is the time at which inflation begins. The fourier decomposition of is 

d 3 p 

where p is the comoving 3-momentum. The modes p (t) are given in terms of Hankel 
functions [11], whose explicit form we shall not need in this paper, and a p , (a' p ) are 
anihilation (creation) operators respectively. In an infinite de Sitter space, the inte- 
gration over 3- momenta p would continue down to p = 0. However, in the inflationary 
Universe there is the inflationary horizon whose comoving size Hj _1 , acts as an infrared 
cutoff on the spectrum. 

Next we split <fi into two parts, one which includes averaging over modes with 
physical momentum k = pexp (—Hit) > Hi (where H, is the physical distance to the 
horizon at the start of inflation), and the other a kind of coarse grained field smoothed 
out over horizon size scales, 



0(x,t) = / ^^e(p-6H i e H ' t ){a p p (t)e ip - x + h.c} 

J (27TJ 

+ X(x,t)-x (3) 

In eq.(3) % is the mean value of the coarse grained field x during inflation, whose 
value we shall discuss later, e is a small parameter 0<e<l, which ensures that the 
'filter' function represented by the step function in eq.(3) rapidly vanishes for comoving 
momenta p »s Hje Hit [11]. By solving the equations of motion for the modes p , in 
the linearized approximation, one may deduce that slowly varying fields x satisfy a 
Langevin type equation, 

X(x, t) = f(x, t) + (^e- 2H ^ V 2 X - ^) (4) 

2 Here we assume that the field responsible for inflation is not (f>, but some other field in another 
sector of the theory. 
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where the noise f(x, t) is given by 

f ( x > = / 7^372 6 ( P " eeHitH i ) e eHit H i 2 { a P P (t)e ip - x + h.c } (5) 
J (2ir) 

Since [11], one can also prove that the noise f(x, t) satisfies (f(x, t)f(x, t)) = 
(Hj 3 /47r 2 )5(t — t), one may derive a Fokker-Planck (F-P) equation from the Langevin 
equation (4). Introducing the normalized probability distribution P(x, Xi t) f° r the field 

X 

/+oo 
dxF( X -x)P(x,X,t) (6) 
-oo 

where F is an arbitrary function, and ( ) f signifies averaging with respect to the noise 
f defined earlier, the corresponding F-P equation takes the form 

The boundary condition on P(x, x,t) at t= (t = t ) , the start of inflation, is taken 
to be P(x,x, t) = S(x — x)- With this initial condition, P(x,x, t) for t > can 
be interpreted as the fraction of the initial comoving volume occupied by the coarse 
grained field x & t time t. 

We are interested in scalar field potentials of the form 

V( X ) = m x 2 M 2 [cos(^) + l] (8) 

which occur in models describing light domain walls [1]. Since in these models the 
mass scale Hj, whilst m x <Mwe suspect that solution for the F.P. eq.(7) is well 
approximated by the solution to the equation with V = 0. To see this, (following the 
derivations given in [11] for general V), we can obtain an approximate solution to the 
F-P equation, with V given in eq. (8) as 

with r = Mt/47r, and where the dispersion A 2 (t) satisfies 

dA 2 (r) 87rm 2 . , a9/ . h 3 

-*r = 3h# cosi(r)A (T) + v < 10 » 

In eq.(10) h = Hj/M, and x(r) is an intermediate variable satisfying 
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. , . 47rm x 2 

x(r) = - * smx(r) (11) 
v ) 3hM 2 v ; v j 

The point is that from the solution to eq.(ll), i.e. 

i / J . .47rrm v 2 . , , , 

x(r) = 2tan- 1 (tan(^) exp (^f ) ) (12) 

we see that since the ratio m x /M <cl [1], x(r) ps x , for all relevant time scales during 
inflation, where x is a constant of integration. This allows one to solve eq.(10) directly 
for the dispersion 

a2/ , 3h 4 M 2 . r 87rm 2 . . . , . 

87T z m x / cos(x ) 3nJvr 

Here A(r) satisfies A(r = 0) = corresponding to our previous boundary condition 
P(x> X? t) = 5(x — x) at t = 0. However, again we see the very small factors of m x 2 /M 2 
in the exponent of eq.(13), so we can expand and keep lowest order terms only. This 
gives 

A 2 M^t + 0(^) (14) 

Hence, for the scalar potentials of relevance to theories describing light domain walls, 
it is a very good approximation to take as the solution to the F-P equation (7) the 
usual Gaussian probability distribution 



^"•"iS^'-T 1 (15) 

centred about x — X- 

Thus for example the correlation of x(x, t) with itself is 

(( X (x, t) - x)(x(x, t) - x)) = + 0((^) 2 ) (16) 

Even if one had considered the correlation function in eq.(16) at two separate (co- 
moving) points x, y the r.h.s. is typically only small at extremely large physical sepa- 
rations | x — y | expHjt ~ Hj _1 expHjt. Because the physical correlation length of x is 
thus typically of order the inflationary horizon, x can be interpreted as a quasi-classical 
field on very much smaller distance scales [11] . 
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3. Evolution of p(x,x) to a non- Gaussian probability distribution. 

As we have indicated above and as is explained in ref. [4], the fluctuations about the 
mean value of x gi ye rise to an inhomogeneous quasi-classical field in the Universe after 
inflation has ended, and the Robertson- Walker (R-W) phase begins. As explained in 
ref. [12], the interpretation of this field as being quasi-classical is clear when thinking 
in terms of physical momenta and length scales, as opposed to comoving quantities. 
The amplitude of the fluctuations in this field, 5x 2 , on physical length scales / in the 
R-W Universe is given by 

^ W = S"4 (17) 

where l c represents an effective ultra-violet cutoff on the physical momentum of the 
modes <p v which contribute to the amplitude in eq.(16). The formula for the dispersion 
o given in eq.(17) above, is clearly related to the t- dependent fluctuations (eq.(16)) 
produced during inflation. The origin of the t-dependence of the r.h.s. of that equation 
ultimately comes from modes with physical momenta k ~ H^exp — Hjt. 

The explicit form of eq.(17) may be understood as follows. After inflation ends 
and the R-W phase begins, new particle horizons appear, whose size are many orders 
of magnitude smaller than that of the 'blown-up' inflationary horizon. In fact the 
presently observed Universe can fit many times over inside the inflationary horizon. 
So we have to reconsider which modes contribute to the fluctuations in x when we 
observe on a physical length scale I which is bounded by the presently observed horizon. 
Obviously this last restriction means that we should only keep modes whose physical 
momentum k satisfies 2n/k < I, so / acts as an infrared cutoff on the spectrum of P . 
At the other extreme, modes whose physical momentum is sufficiently large so that 
2n/k < hU 1 have decaying amplitudes (here H is the Hubble constant during the R-W 
phase ) and do not contribute significantly to the fluctuations. Therefore l c H -1 . By 
contrast, the modes whose wavelengths are greater than the R-W horizon have constant 
amplitudes. All of these properties can be deduced by studying the time evolution of 
the P in R-W space. 

To summarize the above, the emergence of R-W phase places effective cutoffs on 
the spectrum of P . Taking these into account exactly reproduces the amplitude of 
fluctuations as given in eq.(17). 

Before moving on to discuss the behaviour of these fluctuations at later times in the 
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R-W Universe, let us comment on the possible values of the quantity x.(The following 
arguments were originally used in discussing the mean value of the Peccei-Quinn axion 
field [13]). If the field x had appeared after inflation, then microphysical processes 
would tend to smooth the field out over horizon volumes with, in general, different 
values of x in different causally disconnected regions. (Here we are discussing an earlier 
epoch of the R-W phase when the cosine potential can be neglected). In this situation 
it would be natural to take x to be the r.m.s. of all the possible values between — n 
and + 7r namely Xrms = n/V3 [13]. If, on the other hand x is present during inflation, 
(which is the case we consider in this paper), a quite different conclusion holds. The 
point is that a single inflationary patch easily encompasses our present Universe, so 
that the value of x i n this patch is the one we see. But this value is actually arbitrary, 
because whilst the average of x over & U inflationary patches would again be the r.m.s. 
value ir/y/3, this says nothing about the value of x m an y one specific patch. 

This result has important consequences for the formation of light domain walls via 
the fluctuations described previously. When we try to analyze the wall distribution 
in terms of percolation theory, x is related to the percolation probability p[3] (to be 
discussed later). The fact that this is now arbitrary, leads to a richer (and more 
realistic) class of wall structures than had previously been thought possible [9]. 

Associated to the R-W space fluctuations in x there is a corresponding probability 
density P(x,x, l/l c ) given by 

/-\-oo 
dx(x-x) 2 P(x,xJ/Q (18) 
-oo 

which is just the Gaussian distribution of eq.(15) with de Sitter time t replaced by 
Hj _1 ln(/// C ). It is the distribution function in eq.(18) that is important in determining 
the probability of wall formation. In particular wall creation becomes significant if the 
dispersion a is comparable to M [3]. 

If one is to go beyond an estimation of the probability of wall formation, and 
discuss the kind of structures the walls create immediately after formation, statistical 
percolation theory seems to offer some hope. For example in [9], it was shown that 
assuming domain wall distribution can be described by that of percolation clusters in a 
percolating network, hierarchical matter clustering occurs. If matter clustering leads to 
galaxy formation, then percolation theory predicts galaxy clusters of specific richness 
(total number of galaxies) and distribution which fits fairly well with those catalogued 
by Abell [5]. 
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Having said this, one really has to justify applying percolation theory to this prob- 
lem. In 'two-colour' percolation, a lattice site (for example ) is occupied or unoccupied 
with probability p, and 1 — p respectively. The analogue of this lattice in the cosmolog- 
ical case, is the R-W Universe at some time t, being filled with causally disconnected 
regions of size H _1 (t). Each of these regions contains a value of x so that over the 
whole Universe, they are distributed according to the Gaussian P(x, X-, l/l c ) °f eq.(18). 
Clearly P(x, X, l/h), cannot as it stands, give a distribution of x values which are just 
±7T everywhere. In the remaining part of this paper, we will show however, that this 
situation will change if we take into account the effect of the cosine potential on the 
values of x- We will then find that the Gaussian distribution will evolve into one which 
is non-Gaussian, displaying peaks near the two minima x = i 71 " °f the potential. In 
this way we have a situation which is much closer to the one relevant to percolation 
theory. 

To see how P can evolve into non- Gaussian form it is simpler to think of a discrete 

set of distinct random field values X^t), X r (t) 5 which occupy N causally separate 

regions at time t. (N is like the lattice size in percolation models ). The distribution of 
such values will be given by a discrete version of the gaussian P, so long as the values of 
X s (t) for s=l..r remain relatively constant in time. Since each x s (t) satisfies an equation 
of motion like that in (1) with V given by eq.(8) and scale factor a(t) appropriate to the 
R-W Universe, the latter situation will be the case when the friction due to the Hubble 
expansion is greater than the acceleration on each x s due to the cosine potential. That 
is, P will maintain its Gaussian form for t < t , where t is the time at which the x s 

begin to roll down their potential. We shall refer to the values x^t), X r (t) for t 

< to as Xo, ■■■Xo- K n i5 ••• n r are the number of causal volumes occupied by the values 
Xo, ■■■Xo then the probability of finding the value xl f° r s — 1> - r is n s/N . 

Now let us consider what happens when t > to- We have to consider solutions of 
eq.(l), with the boundary conditions x s (t = t ) = Xo an d X s (t = t ) = 0, s=l,..r . The 
solutions x^t)? --X r (t) are trajectories with initial points xl, ---Xo? which in general, will 
cross each other at some time t = t c > t (whose value we shall discuss later) even 
if the initial values are all distinct. This is because the equation of motion for each 
X s form a non-autonomous set of differential equations which depend explicitly on time 
t. However, for to < t < t c , the trajectories do not cross and we can approximate 
them as solutions to eq.(l) in the overdamped regime where the term involving second 
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derivatives is dropped. The resulting equation is then easily seen to be autonomous, 
so this procedure is consistent. 

The point is that since the trajectories x s (t) are unique for t < t < t c , the proba- 
bility of finding the value x s (t) will be the same as for t < t i.e. for the value xl- This 
fact will allow us to determine the functional form of P t on x(Xo,t), where we now 
revert back to continuum notation, and P t denotes the distribution of the t-dependent 
field x(xojt) for t < t < t c . That is 



Pt(x(Xo, t), Xo, v) = P(Xo, Xo, o") 
Solving eq.(19) for P t , we find 



(19) 



P -J 2ir ^( 2M V(x(X0,t)) xo_ 2 
1 " \l .2 ex P i 2M 2M Jj 



(20) 



where 



f(x(Xo,t)) = Xo (21) 

i.e. f is the inverse function that maps the solution x(Xo,t) — > xo- Here we see the 
importance of requiring the trajectories be unique, otherwise the inverse mapping f 
would not exist for all initial values Xo- 

Finally, we have to take care of the correct normalization of P t . Eq.(19) was moti- 
vated from the discrete approach, whereas in passing to the corresponding continuum 
expressions we have to worry about Jacobian factors in going from P(xo) Pt(x(t))- 
The latter involves changing variables from xo — f(x(t)) to x(t) 

roo roc (Jf 

y_ oo dxoP(xo,Xo,^) = l= /^dxW^yPtWt),^,^) (22) 
Hence, the correctly normalized continuum distribution is 

Pt(x(t),*o,<7) = exp(-^[f(x(t)) - Xo] 2 ) (23) 

To see what P t looks like, we need to solve eq.(l) in the overdamped approximation, 
and then compute the inverse function f. The solution for x(t) satisfying the boundary 
conditions discussed earlier, is 
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X (t) = 2Mtan- 1 (tan(||)e A ^')) (24) 
m 2 1 1 

A{t - t ° )= ^m-wM > (25) 

It is then straightforward to obtain the function f, which gives the following P t 



Pt(x(t),Xo,0") = 



27T e -A(t,t ) 




a 2 ( cos 2 (|g) + sin 2 (U ) e ~ 2A (^y 



x exp _(^[tan- 1 (tan(^)e- A ^)-|^] 2 ) (26) 

Let us study the distribution P t for the simple case when xo — i.e. when the original 
Gaussian is symmetric about Xo — 0, as this will allow us to see the departure from 
Gaussian behaviour. We look first for the stationary points of P t , which satisfy 

tan(gl) = tan[ ^ a 2 sin(#) ] e A (27) 

v 2M ; L 4M 2 M y J 

Obvious solutions to eq.(27) are x = mod 2-7T, but there are other solutions that 
depend on A but which are difficult to find analytically. However, let us concentrate 
on the central stationary point at x — 0- One can show that 

/d 2 P e A -e- A 1 , 

sl s n ( = Slgn( ~2yp~ - ^ } (28) 

Thus, at t = to, A = and the point xo — is clearly a maximum, as it should be since 
P is a Gaussian centred at zero. As A increases with time, this situation continues 
until a critical value t cr i t is reached , which is defined by 



e A cri t _ e -A crit = ^ 



2M 



a 



When A > A cr i t , x — changes from a maximum to a minimum. Now, on the other 
hand, P t is a bounded function as x(t) — > ±oo hence if there is a minimum at x = 
there must exist other maxima (at least 2 ) at some non-zero values of x- We expect 
that similar conclusion holds for distributions with x ^ 0, which although difficult to 
prove analytically, is verified by plots of P t (see figs.3 ). 



11 



Before discussing these numerical plots, we should return to the issue concerning 
the existence of the crossing time t c beyond which trajectories x(Xo,t) with different 
initial values Xo> cross. In fig.l, numerical solutions to the full equations of motion 
are plotted as solid curves for various initial values Xo- It is clear from these, that t c 
exists, and is approximately given as m x t c ~ 3. Note we have not considered initial 
conditions beyond Xo/M = ±7r because P t is exponentially suppressed in this regime. 
By contrast, the dashed curves in fig. 2 are the analytic solutions given in eqs. (24,25). 
We see that for m x t < 0(1), the latter are a good approximation to the numerical 
solutions of the complete equations of motion. Thus, the expression for P t of eq.(26) 
is only correct for times up to maximum t max , with m x t max ~ 0(1). 

In figs.3a-3d we have plotted the distribution Po = Vt/ ^2ti /a 2 for values of 
(x,<7/7rM) corresponding to (0,0.5), (0,0.8), (7r/4 , 0.5) and (7r/4, 0.8) respectively. 
The plots are for time t \^2m x ~ 1 . The values of a/nM are chosen so as to display 
the double maxima in P t , described earlier. Such values are very reasonable, because 
as was shown in ref. [3], the probability of domain walls forming in these models can 
be estimated as the ratio cr/irM. The dashed curves in figs. 3 are those of the original 
Gaussian distribution at t = t . It is clear from these plots that the effect of the cosine 
potential is to focus probability into the minima of the potential at x = i 71 ") an d we 
can see roughly that this probability is equal to the area under the dashed curve to 
the right or left of the point x — Xo- These two areas play the role of the probability 
p and (1 — p) in percolation theory. 

In conclusion, we have shown how the probability distribution P(x, x) generated 
during inflation evolves from Gaussian to non-Gaussian behaviour during the R-W 
phase. This was a consequence of the cosine potential of x focussing probability into 
the two minima at x — ^ zn - The amount of focussing in each minima depends on the 
value of x- These results justify the percolative approach to domain wall formation 
in this model, with the important consequence that the percolation probability p can 
take any value < p < 1. 
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